Sobolev spaces and Lagrange interpolation 



B. Bojarski* 

Institute of Mathematics, Polish Academy of Sciences 
00-956 Warszawa, Poland 
E-mail: b.bojarski@impan.pl 

January 24, 2012 



In this short paper the discussion of the pointwise characterization of func- 
tions / in the Sobolev space W m ' p (W rl ) given in the recent paper [6] is supple- 
mented in §1 by a direct, essentially geometric, proof of the novel inequality (for 
m > 1), appearing in [6 apparently for the first time, and involving the use of 
the m-th difference of the function /. 

Moreover in §2 some additional comments to the text in [6] are given and a 
natural class of Sobolev spaces in domains G in 1" is defined. §3 contains some 
final remarks. 



Let us recall [15] , [19] that for an arbitrary integer I > and a real or complex 
valued function / on R" the expression 

A^):=E(- 1 )^(0 /(x + ^ ) = ( " 1)i i ( " 1) '( / ) /(a;+j/l) (L1) 

is called the Z-th difference of the function / at the point x 6 M. n with step h, 
h G K", h ± 0. We set also A° h f(x) := f(x) and A l f(x) := 0. As is classi- 
cally known, for y = x + Ih, h = ^j^, (11.11) has a beautiful interpretation as 
the difference or error in approximating the function f{y) by its interpolating 
polynomial evaluated at y 

A l f(x; y) = A{f(x) = f(y) - C(y; f;x ,..., x^) = (-l) l A l f(x, y) 

{x =x), (1.2) 

where C(y; f;xo, . . . , xi-\) = 2j=o f{ x i)^j(V-> x o, ■ ■ ■ is the Lagrange in- 

terpolating polynomial for the function / and the equidistant colinear nodes 
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Xi = xq + ih, i = 0, ...,/ — 1, [TS]. Here £j(y, xq, • ■ ■ , xj_i) stand for the funda- 
mental Lagrange polynomials in y. A l (x,y) is the notation used in [BJ. Let us 
remark that all points Xj r = x + jh, j = 0, . . . , I — 1, X; = t/, are situated on the 
afhne line 7?. in R™, joining x and y, which can be identified with the real line 
R 1 , making all the algebraic operations inherent in (jl.2p meaningful. 

In the sequel we shall use for (jl.2l) the term Lagrange interpolation remain- 
der, or just Lagrange remainder, of order I at the point x evaluated at y, in 
analogy with the term Taylor-Whitney remainder centered at x, 

R l - 1 f{x;y):=f{y)-Tlr 1 f{y), (1.3) 

now common it mathematical literature. In [BJ A'/(x, y) was also called an Z-th 
finite difference remainder of the function / at x evaluated at y. 

For functions / 6 W m ' p (M. n ), p > 1, the fundamental novel inequality re- 
ferred to above reads as 

|A™/(x;y)| < \x-y\ m [af(x)+af(y)} (1.4) 

for some a f € L p (R n ). 

We skip here over the somewhat delicate point that the Sobolev functions in 
general do not have pointwise values and the left hand side of inequality (II .4[) 
is meaningful only up to subsets of measure zero. The right hand side may be 
infinite on a non-empty set of measure zero. 

The functional coefficients a™(x) in (|1.4p are not uniquely defined. They 
are collectively called mean maximal m-gradients of the function / and play 
the role of a variable Lipschitz coefficient of /. Roughly speaking, they all can 
be majorized by the local maximal function of the generalized Sobolev gradient 
|V m /| of / e W m ' p (R n ) as will be also seen from the constructive proof of (|1.4D 
sketched below. 

Our proof is organized in a series of lemmata. 

Lemma 1. For f £ W ,P (R"), 1 < p < oo, the following inequality holds 

\f(x)-f(y)\<\x~y\(a s f (x)+a 5 f (y)), x, y e R n (1.5) 

for some a s f 6 Lf oc (K"), 6 = \z - y\. 
Proof. For arbitrary x, y G R" we have 

f(x)-f(y)=[(Vf(x + ht),h)dt, h = y-x, (1.6) 
Jo 

hence ^ 

\f(x)- f(y)\<\x-y\ f \Vf\(x + ht)dt. (1.7) 
Jo 

Let B(x,r) be the ball of radius r centered at x, and E r (x, y) the spherical 
segment 

S r = B(i,r)n%,r), r = |x-y|. 
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For an arbitrary z G S r 

\f(x) - f(y)\ < \f(x) - f(z)\ + \f(z) - f(y)\. (1.8) 
Since \x — z\ < \x — y\, \y — z\ < \x — y\, averaging (|1.8[) over z G S r we get 

\f(x) - f(y)\ < l\f{z) - f(x)\ da z + -f\f(z) - f(y)\ da z 



j \f(z)-f(x)\da z + j \f(z)-f(y)\da z ), (1.9) 

B(x,r) -B(j/,r) 

where the notation \G\ for a subset G in R" is used for the volume of G, 
\B(x,r)\ = \B(y,r)\. 

By elementary geometry the ratio ^^f'j"^ is a constant depending only on n, 

^'j^ = C(n), and, as is well known, the average j- g ^ x r \ \f(z) — f(x)\ da z is 
estimated by the local Hardy-Littlewood maximal function at x of the gradient 
|V/|, M 5 (\Vf\)(x), gg, g3]. Here the inequality (TO]) is used. Thus in (TT3]) the 
function a* (a;) is controlled by M 5 {\ V/|) (x): in fact a£ (a) < C(n)M 5 (\7f\)(x). 

□ 

The proof above, without changes, works for vector valued functions. This 
proof should be compared with the proof of the basic pointwise inequality (1) in 
our paper with P. Hajlasz from 1993 [7]. Notice that it does not refer to Riesz 
potentials and Hedberg lemma as in [7]. 

When combined with Reshetnyak's trick [36] used in [7], it can be used 
to deduce, in a direct way, a new and simplified proof of the basic pointwise 
inequalities in [7J [S] . 

Let / 6 W k ' p (R n ). For / = 0, 1, . . . , k - 1, h G M", consider the functions 



9 l h {x)= [ ■■■ [ V l f(x + y2t i h)(h,...,h)dt 1 ...dt l , (1.10) 
Jo Jo K i=i ' 



where V / is the Z-th gradient of / considered as an ?-poly linear form on M. n . 
Lemma 2. g l h (x) as a function of x G K n is in the class W k ~ Lp (R n ). 
Proof. Obvious. □ 
In particular, <?°(x) = f(x), 

9h(x)= [ (Vf(x + th,h)dt= [ V ^-(x + ht)hidt, etc. 
Jo Jo trf dx % 
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Lemma 3. The function g l h (x) for 1 = 0,..., k—1 has the integral representatic 



9h 



3=0 



[ ••• f V l f(x + Y^Uh)(h,...,h)dt 1 

JO JO V i=1 ' v ' 

t times 



dti. (l.ii) 



Proof. This is the well known formula of finite difference calculus [T^l [TH] . For 
I = k it is also used in the paper [TT] of R. Borghol. □ 

Lemma 4. For I = k — 1 we /icwe i/ie formula 

9t\x) ~ 9h~H* + h) = E(-l)'- 1 (* Zl) [K x + V- W ~ f( x + lh )] 

i=i ^ ' 



3=0 

in the notation of [6 . 
Proof. By Lemma [3] 



k 



Hence for ft, = ^j- 



A k h f(x)=g k h - 1 (x)-g k h - 1 (x + h) 



2=1 

fc-1 



E(-l)' -1 ( ! _ ! ) + V ~ W - /(* + "»)] 



Ec-iy'r . 1 )/(^+^)+E(- 1 ) i (z: 1 1 )/( a: +^ 

j=o \ J / l=1 V / 



= /(*) + (-l) fe /(* + fc/0 + E(-!) J [(* j l ) + Q _ + 

= (-l) fc A^/(z,y)=A fe /(*,y). □ (1.14) 

Formulas (|1 . lpp — (|1 . 14D are examples of a series of formulas of finite differ- 
ences [T3], Q35] which connect the operations of vector differential operators V l f 
with finite difference operators Ixif for various values of the parameters I, k 
and h. They allow to reduce the estimates of higher order difference remainders 
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of functions to lower order remainders of higher order gradients of these func- 
tions. They are analogues to the operations in Taylor- Whitney's algebras which 
played an important role in the development and applications of pointwise in- 
equalities in [2] , [5] , [BJ , [5] and in the Whitney-Glaeser-Malgrange theory [2U] , 
[25] . [5T] . [35] . of smooth functions on arbitrary closed subsets of R". Deep 
and important papers of Glaeser [2T], [32], [33] in this theory seem to be so far 
waiting for better understanding and exploitment. 

Now, combining Lemmata [T]-[5] we obtain the estimate (|1.4j) . with a/(x) 
controlled by the maximal function of the vector gradient V k f as required. 

For convenience we summarize our discussion in the following 

Proposition 5. Let f 6 W m,p (M. n ), 1 < p < oo. Then there exists a function 
cif £ L p (R n ) such that the inequality (|1.4I) holds for almost all points x, y € M. n . 
The function cif is majorized a.e. by the local maximal function of the generalized 
Sobolev gradient |V m /| € L p (R n )\. 

As formulated above Proposition [5] is the "necessary part" of the main the- 
orem in [5J for the case of Lagrange interpolation remainders. The proof of 
the "sufficiency part" of the theorem is left unchanged and proceeds along the 
argument sketched in [5J. 

The sketch of the proof of Proposition[5]presented above is in the convention 
that we work in the class of smooth functions where all the operations involved 
are classically meaningful. The clue of the story is that the constants appearing 
in the estimates depend on the parameters n and p only. This fundamental fact 
comes up again in §2 below where the pointwise inequality appears as "stable" 
under smoothing by convolution. 



One of the advantages of the Lagrange interpolation calculus and the Lagrange 
remainders over the Taylor- Whitney remainders is that they interact very well 
with convolutions. This is immediately seen for the remainders A 1 f(x; y) = 



Indeed, for a normalized mollifier ip £ (x), ip e > 0, J <p s (x) dx = 1, we have 
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f{y)-m- 




and 




hence 




a(x - rj)(p e (rj) d-q + a(y - r])<p e (r)) drj 



or 



fe(y) - fe(x)\ < \x-y\(a e (x) + a £ (y)). 



(2.1) 
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This elementary though basic fact, already formulated in our paper [7] and 
even earlier and repeated later in many seminar talks, holds for the higher 
order remainders A m /(x; y) as well. Indeed, e.g. for m = 2 we have 

f e (x) - 2f e (2p) + f £ (y) = J [fix - v) ~ 2/(4* - V ) + f(y - n)]<p e (rj) dr, 

< \x-y\ 2 J [a f (x - rj) + a f (y - rj)]ip e {rj) dr] = \x - y\ 2 (a f , £ (x) + 0/, e (y)) 

and the same calculation works for m > 2. 

Pointwise inequalities for the Lagrange remainders A m f(x, y) appeared also 
in the recent papers of H. Triebel and his school [55] , [IS] . When the paper [5] 
was written the papers of Triebel [45] and Haroske-Triebel [26] were unknown 
to the author. Geometrically for m > 1 they differ from ours by introducing the 
intermediate nodes Xi = x + ih, i = 1, . . . , m — 1, in the right hand side of (|1.4p . 

For arbitrary 0<p<oo,s>0 and m € N with s < m in [45] is introduced 
the class Lp' m (R™) of all / € L p (W l ) for which there exists a nonnegative 
function g £ L p (R n ) such that for all h G R™, < \h\ < 1 the inequality 

m 

|A™/(x)| < \h\ s 9{x + lh) a.e. in W 1 (2.2) 

1=0 

holds. 

With the norm 

||/lkm,p = H/HmR") + illf HfllipfK"), 

where the infimum is taken over all g admissible in (|2.2[) . the space Lp' m (R n ) is 
a quasi-Banach space [45] . 

For s = m, p > 1, the inequality (I2.2[) is stronger than (|1.4I) . i.e. (|1.4j) implies 
(|2.2p . Now it is clear that one of immediate conclusions from the main theorem 
in [6] is that for these values of s and p (|1.4[) is equivalent to (|2.2[) . Thus we 
conclude that the intermediate nodes in (|2.2j) can be discarded. 

For a while, let us introduce the notation: W m ' p (R") is the class of all 
functions / G L P (R") for which there exists an a/ £ L p (R") such that inequality 
CL3} holds a.e. in R™. 

In the general context of (|1.4I) and (12 .2[) we have the following 

Proposition 6. For f £ W m ' p (R ra ) the mollified function f £ — f * ip £ is in the 
class W m ' p (R") and 

\A m f £ (x;y)\ < \x- y r(afjx)+afjy)) (2.3) 

with a™ e (x) — a™ * ip E . 

Moreover, by the known properties of convolutions the L p -norms of the mean 
maximal gradients aj" e are uniformly controlled by the L p -norms of mean max- 
imal gradients of f: 

ll^ellMR") < WW*") ( 2 - 4 ) 

for all e > 0. 
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In particular, we conclude that smooth functions are dense in W m ' p (R n ). 
For the same values of the parameters s,m,p as for Lp' m (R") we can also 
consider the class Lp' m (R") defined by the pointwise inequality 

\A m f(x;y)\ < \x-y\ s (g(x)+g(y)), y = x + mh, (2.5) 

for some g £ L p (R n ). 

In [IS] the quasi-Banach spaces Lp' m (R n ) are used to identify some Besov 
spaces Bp 6 q (R n ), O<0<l,O<g<oo, as real interpolation spaces ([35], 
the main theorem). For s = m in [B], as well as in the paper |45] . the spaces 
W" l ' p (R™) are identified with classical Sobolev spaces W m ^(W l ). 

The natural interesting question is to characterize the spaces Lp' m (R n ), for 
s not integer, as some Sobolev-Besov type spaces. 

Proposition [6] and its proof are also valid for the class Lp' m (R"). Thus 

smooth functions are also dense in Lp' m (R n ). 

The described characterization of Sobolev spaces W m,p (M. n ), p > 1, by point- 
wise inequalities (jl.4[) obviously holds for open subdomains G C R" as well, if 
they have sufficiently regular boundary, e.g. for extension domains [38 . How- 
ever this characterization is definitely not true for arbitrary subdomains. In 
this context it seems legitimate to introduce the (maximal) class of subdomains 
which admit global pointwise characterization. 

Definition 7. An open subdomain G C M™ is called a natural Sobolev (p,s), 
1 < p < oo, s > 0, domain if the pointwise inequality 

\A m f( x> y)\<\x-y\ s (a(x) + a (y)), a < m, (2.6) 

for all pairs of points x, y £ G such that the segment [x, y] C G, defines a 
Sobolev type Banach space W S ' P (G). 

For s < m, 0<p<ooa related class of spaces (K") m has been introduced 
by H. Triebel in [45] and identified with subspaces of Besov spaces Bp ^(K™) 
for G = R". 

It is natural to ask in particular in what sense our Sobolev-Besov type spaces 
coincide with Triebel-Besov spaces in [45]. More generally we can ask in what 
sense and for what subdomains the spaces W S,P (G) coincide with Besov-Sobolcv 
type spaces for non-integer s. 

3 

The pointwise inequality (|1.4[) actually can be used to characterize rather the 
homogeneous Sobolev spaces W m,p (G) for a subdomain G C R™ with the semi- 
norm || V m / 1| lp(g)- The classical inhomogeneous Sobolev spaces arise then as 
subspaces of L P (G). For the model case G = M™ various delicate phenomena 
are related to the asymptotic behavior of functions in W m ' p (M. n ) for |a;| — > oo 
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which seem to be so far only partially understood (see [3] and numerous other 
references. See also [1] for the case n — 1, i.e. on the line K 1 ). 

In our presentation here the pointwise inequality (11.41) as well as some- 
what more sophisticated inequality for Taylor- Whitney remainders R™^ 1 f(x\ y) 
should be considered as elementary, though fundamental, facts appearing at the 
first steps of any discussion of Sobolev spaces. 

What seems to be still lacking in this elementary discussion of Sobolev space 
theory, is the deeper, geometric and analytical, understanding of the trace (pro- 
jection) operator W m,p (W l ) — > W s ' p (R k ) for the corresponding values of the 
parameters s,p. As is classically known, this question led to the introduction of 
Sobolev fractional spaces, Vy s < p (R™), s — real, H s - p (R n ) and Besov spaces ([3J, 
[42], [43] and many references therein). 

It seems also that special attention should be directed to u les schemas 
d 'interpolation" of G. Glaeser [3U], [53], [37]; [31], and their role in the the- 
ory of Sobolev spaces, see also [37], [38] , 

As already remarked in [BJ, and earlier even in [3], [7J, [5J, the pointwise in- 
equality (|1.4[) above, together with the related inequality for the Taylor-Whitney 
remainder R m ~ 1 f(x;y) (precisely, inequality (1.2) in [6j) may serve as natural 
starting points and effective tools in the study of fundamental structural prop- 
erties of Sobolev functions. Let us briefly recall some of them without going 
here into details (postponed to the activity and exposition plan foreseen in the 
last lines of [BJ). 

a) Lusin's approximation of Sobolev functions, i.e. interpolation by smooth 
functions on closed subsets, up to open complements of arbitrary small 
measure. 

b) Stability under the convolution with compactly supported C°° kernels; 
density of subspaces of smooth, C°° functions. 

c) S. M. Nikolskii's [3J fundamental theorems describing the characteriza- 
tion of Sobolev functions by their behavior on typical (almost all, in some 
natural sense) hyperplanes of positive codimensions less than n. The in- 
equalities (|1.4p reduce the characterization of Sobolev functions to their 
behavior on affine segments in their domain of definition. In particular 
the functions in W m ' p (R n ) for p > 1 are Holder continuous on a.e. hy- 
perplanes M fc C R™, k < p, the only global condition binding the variable 
Lipschitz coefficients on hyperplanes is the Fubini theorem for Lebesgue 
spaces L p (R n ) and the factorizations R™ - R n ~ k x R k . 

d) Characterization of compact subsets of W m ' p (R n ) by some conditions on 
the corresponding variable Lipschitz coefficients. 

e) Differentiability properties of Sobolev functions, Calderon differentiability 
theorems, approximate and Peano differentiability, [3J, [7J, [5], [T3J, [T4"] . 
ES], [33], [35], 02. 
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f) Extension of Hajiasz-Sobolev imbedding theorems of Sobolev spaces into 
higher exponent Lebesgue spaces: W l ' p (G) C L q {G), q > p, I > 1, for 
suitable values of the parameters l,p, q and dimG = n, as in the classical 
Sobolev theory, modeled on Hajlasz's proof for measure metric spaces, 

g) Extension of the classical Hermitc interpolation formulas [TS], [TO] to the 
general multidimensional context of "multiple" nodes as a theory inter- 
mediate between the Taylor- Whitney and Lagrange ("simple" nodes) in- 
terpolation theory ([20 ] . [25 ] . [3T]). 

h) If instead of the colinear equidistant nodes other configurations of inter- 
polation nodes are used, more complicated algebraic and geometric phe- 
nomena occur, e.g. in the case of colinear not equidistant simple nodes the 
classical divided difference calculus comes up fTJ] , [TB] . New interpolation 
concepts appear in the Glaeser papers [20], [22], and, more recent, [23] , 
|31) . Combining these ideas with the Sobolev's averaging procedure, [39] . 
[40) , [41) , defines an apparently new interesting research direction. 

i) Also the natural inclusions, e.g. W m+1 ' p (R n ) C W m ' p (M. n ) when inter- 
preted in terms of pointwise inequalities (|1.4I) suggest direct implications 
between inequalities (|1.4j) for various admissible values of the parameters 
(m,p). These lead to interesting and non-trivial arguments of geometric 
and analytic character. Probably the first beautiful example of this type 
of argument was given by Y. Zhou [53] . 

Let us remark at last that the "pointwise" approach to the Sobolev space 
theory seems to bring out more clearly than usually presented in the literature 
(e.g. [30], [32], [33J, [TS], [I3J) close and natural connections between the general 
concepts of functional spaces and approximation theory on the real line R 1 and 
in n-dimensional, n > 1, euclidean spaces. 
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